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Abstract 

Inspired by the recent developments in the fields of quantum distributed computing, 
quantum systems are analyzed as networks of quantum nodes to reduce the complexity of 
the analysis. This gives rise to the distributed quantum consensus algorithms. Focus of 
this paper is on optimizing the convergence rate of the continuous time quantum consensus 
algorithm over a quantum network with N qudits. It is shown that the optimal convergence 
rate is independent of the value of d in qudits. First by classifying the induced graphs as 
the Schreier graphs, they are categorized in terms of the partitions of integer N. Then 
establishing the intertwining relation between one level dominant partitions in the Hasse 
Diagram of integer N, it is proved that the spectrum of the induced graph corresponding 
to the dominant partition is included in that of the less dominant partition. Based on this 
result, the proof of the Aldous’ conjecture is extended to all possible induced graphs and the 
original optimization problem is reduced to optimizing spectral gap of the smallest induced 
graph. By providing the analytical solution to semidefinite programming formulation of the 
obtained problem, closed-form expressions for the optimal results are provided for a wide 
range of topologies. 

Index terms — Quantum Networks, Distributed Consensus, Aldous’ Conjecture, Optimal 
Convergence Rate 

1 Introduction 

Coordinated control and consensus are the essential factors in coupled dynamical systems that 
are modeled as networks of autonomous agents. Examples of such systems in nature are flocks of 
birds, school of fish, neurons within the nervous system [18]. Other manmade examples include 
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power grids and social networks [Uj. Within the fields of distributed control and optimization on 
networks, there have been remarkable efforts in modeling and controlling cooperative collective 
behaviors in networks of autonomous agents. One of the basic behaviors of autonomous agents 
is reaching consensus, where agents in a network achieve a common state using only local 
communication with each other [231II51130]. 

Due to recent progress in the fields of quantum information science and quantum distributed 
computing [SEIE], quantum systems are analyzed as networks of Quantum nodes. This is 
because building and maintaining a centralized and rather big quantum computer with many 
qubits can be excessively difficult and expensive. As an alternative solution, such machine can 
be replaced by a network of smaller quantum computers to carry out computation on a scale 
that would yield practical benefits. 

In |21l [20l |22j authors extended the consensus problem to the quantum domain by reinter¬ 
preting it as a symmetrization problem. They have addressed this problem by a switching-type 
dynamics based on convex combinations of actions of a finite group. Furthermore, they derive 
the general conditions for convergence and show that convergence is guaranteed provided that 
some mild assumptions are held. Using the results on convergence, they prove that it ensures 
asymptotic convergence as well. 

Authors in [271 l2Hl [29] propose a new approach based on the induced graphs of the quantum 
interaction graph for relating the quantum consensus over the A^-qubit network to the classical 
consensus dynamics. They have shown how to carry out convergence speed optimization of the 
equivalent classical consensus via convex programming. Furthermore, they establish necessary 
and sufficient conditions for exponential and asymptotic quantum consensus, respectively, for 
switching quantum interaction graphs. 

In the present paper, we optimize the convergence rate of the quantum consensus over a 
quantum network with N qudits. The main motivation for extending this problem from qubits 
to qudits is to show that the optimal convergence rate is independent of the value of d in qudits. 
By expanding the density matrix in terms of the generalized Gell-Mann matrices, we have shown 
that the induced graphs are the Schreier graphs where in the special case of interchange process 
it reduces to the Cayley graph. Then using the Young Tabloids, we categorize the induced 
graphs obtained from all possible partitions of the integer N. By establishing the intertwining 
relation between one level dominant partitions in the Basse Diagram of integer N, we have 
shown that the eigenvalues of the induced graph corresponding to the dominant partition is 
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included in the eigenvalues of the less dominant partition. Based on this result, we extend the 
proof of the Aldous’ conjecture [Hj to all possible induced graphs and we show that the problem 
of optimizing the convergence rate of the quantum consensus reduces to optimizing the second 
smallest eigenvalue of the Laplacian of the induced graph corresponding to partition {N — 1,1). 
Based on the extension of Aldous’ conjecture, it is shown that if one of the induced graphs 
serves as the underlying graph for another quantum network then the spectral gap and thus the 
convergence rate for induced graphs obtained from the new network is same as those of the old 
one. In this way, a numerous number of weighted Laplacian matrices with the same spectral gap 
can be obtained. In the final stage, we have analytically solved the semidefinite programming 
formulation of the reduced optimization problem for a wide range of topologies and provided 
closed-form expressions for the optimal convergence rate and the optimal weights. 

The rest of the paper is organized as follows. Section presents some preliminaries includ¬ 
ing relevant concepts in graph theory. Young tabloids, Cayley and Schreier Coset graphs. The 
continuous time consensus algorithm and the semidefinite programming formulation of its op¬ 
timization are presented in Section Section describes optimization of the continuous time 
quantum consensus problem and how it can be transformed into optimization of a classical con¬ 
tinuous time consensus problem. In Section analytical optimization of the continuous time 
consensus problem and closed-form expressions for the optimal results for a range of topologies 
have been presented. 

2 Preliminaries 

In this section, we present the fundamental concepts from graph theory, Young tabloids, Basse 
diagrams, Cayley and Schreier coset graphs. 

2.1 Graph Theory 

A graph is defined as ^ = {V,T} with V = {1,..., N} as the set of vertices and £ as the set of 
edges. Each edge {i,j} G £ is an unordered pair of distinct vertices. If no direction is assigned 
to the edges, then the graph is called an undirected graph. Throughout this paper, we consider 
undirected simple graphs with no self-loops and at most one edge between any two different 
vertices. The set of all neighbors of a vertex i is defined as A/) = {j G V : {i,j} G T}. A 
weighted graph is a graph where a weight is associated with every edge according to proper 
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map VF : (S’ — )> M, such that if {i,j} G £, then W{{i,j}) = Wij] otherwise W{{i,j}) = 0. The 
edge structure of the weighted graph Q is described through its adjacency matrix (Ag). The 
adjacency matrix Ag is a N x N matrix with {i,j}-th entry (Ag{i,j)) defined as below 


Ag{i,j) 


Wij if {i,j}e£ 

< 

0 Otherwise 


i.e., the {i,j) — th entry of Ag is 1 if vertex j is a neighbor of vertex i. If the graph Q has no 
self-loops Ag{i,i) = 0, i.e., the diagonal elements of the adjacency matrix are all equal to zero. 
For undirected graphs the adjacency matrix is symmetric, i.e., Ag is symmetric. The degree of 
a vertex i is the sum of the weights on the edges connected to vertex i, i.e. 

N 

di = y^^wjj. 

i=i 


The degree matrix Dg of G is the N x N diagonal matrix where its i-th diagonal element is 
equal to the degree of vertex i and all non-diagonal elements are equal to zero. A graph is called 
connected if there is a path between any two vertices in the graph. A graph is called a regular 
graph if all the vertices have the same number of neighbors. The Laplacian matrix of graph G 
is defined as below, 


Lg{i,j) 


Dg{i,i) if i = j 

< 

-Ag{i,j) if i^j 


This definition of the Laplacian matrix can be expressed in matrix form as Lg = Dg — Ag , where 
Dg and Ag are the degree and the adjacency matrices of the graph G- The Laplacian matrix 
of an undirected graph is a symmetric matrix. The eigenvalues of the Laplacian matrix (Lg) 
are all nonnegative. Defining 1 and 0 as vectors of length iV with all elements equal to one and 
zero, respectively, hence for the Laplacian matrix we have Lg x 1 = 0. In undirected graphs, 
the associated Laplacian is a positive semidefinite matrix and its eigenvalues can be arranged in 
non-decreasing order as follows. 


0 = Ai(Lg) < A2(Lg) < • • • < AAr(Lg) 

The second smallest eigenvalue A 2 (Lg) is known as the algebraic connectivity and reflects the 
degree of connectivity of the graph m- First introduced in m, this eigenvalue is named 
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algebraic connectivity due to its importance in connectivity properties of the graph. Since then 
the algebraic connectivity has found applications in the analysis of numerous problems including 
combinatorial optimization problems such as the maximum cut problem, certain flowing process 
and the traveling salesman problem [T]. The algebraic connectivity can be used to define the 
spectral gap. The spectral gap gives insight into important properties of the graph such as the 
mixing time of random walks [25]. In some cases, the term spectral gap is directly used to refer 
to \ 2 {Lg)- A necessary and sufficient condition for the algebraic connectivity to be nonzero is 
that the graph Q is connected [Tj. If the algebraic connectivity of the graph Q is nonzero then 
Lg is an irreducible matrix i.e. it is not similar to a block upper triangular matrix with two 
blocks via a permutation |14) . The largest eigenvalue X]sf(Lg) of the Laplacian matrix is known 
as the Laplacian spectral radius of Q. 

2.2 Symmetric Group, Young Tabloids, Young snbgroup Hasse Diagrams 

The set of all bijections 11 : {I,-- - ,N} —>■ {I,-- - , N} with composition of maps forms a finite 
group of order A^!, called symmetric or permutation group denoted hy Sn- A standard notation 
for the permutation that sends i —)• n(i) is: 

/ 1 2 3 ... Y X 

Vn(i) n( 2 ) n(3)...n(Y)y 

A r-cycle is a permutation of the form n(/j) = Zj+i for i = 1,... ,r — 1 and n(/r) = h where 
h, - ■ ■ ,lr £ {!) • • ■ ) Y} are distinct from each other and n(r) = z if i not among the Ij. The 
standard notation for this cycle is hr " ) h)- A transposition is a cycle of length 2, and an 

elementary transposition is a transposition of the form (z, z +1). Every permutation 11 G Sjsr can 
be written as a product of disjoint cycles and cycles can be written as a product of elementary 
transpositions. 

A Positive integer Y can be partitioned into a group of positive integers n = (ni, n 2 ,..., nx) 
where their summation is equal to Y and they are sorted in non-increasing order, i.e. ni > tz 2 > 
... > UK- n is referred to as a partition of Y and it is denoted by re h Y. A Young diagram is a fi¬ 
nite set of boxes arranged in left-justified rows with non-increasing lengths. In the corresponding 
Young diagram of the partition re = (rei, re 2 ,..., nx), there are re* boxes in the z-th row of the dia¬ 
gram. As an example, the possible partitions for number 4 are (4), (3,1), (2, 2), (2,1,1), (1,1,1,1) 
and their corresponding Young diagrams are depicted in figure It is obvious that there is a 
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(4) (3,1) (2,2) (2,1,1) (1,1,1,1) 

Figure 1; Young diagram for all possible partitions of 4. 


1 

2 

3 

4 


Figure 2: Young tabloid for partition (2,2). 

one-to-one correspondence between partitions and the Young diagrams. For a given partition 
n h A^, a Young tableau of n-shape is obtained by filling in the boxes of the corresponding 
Young diagram of partition n with integers from 1 to N. For a given partition if the integers 
in rows and columns are ordered in increasing order, then the Young tableau is referred to as 
standard Young tableau. In total for a given Young diagram, there are N\ Young tableaux. The 
definitions above for Young tableau and diagram are adapted from |26j . 

A Young tabloid is an equivalence class of Young tableau under the relation that two tableau 
are equivalent if each row contains the same elements. The notation used for the Young tabloid 
is similar to the Young tableau but without vertical bars separating the entries within each row. 
For a given partition n, the number of Young Tabloids is equal to zz = N\/{ni\ ■ n 2 \ ■ ■ ■ uk^.). 

Young tabloids of a given partition n can be characterized in simple way by so called Ya- 
manouchi symbols. For any Young tabloid of n-shape we define a Yamanouchi symbol as a row 
of N numbers (ri,r 2 ,--- ,rjv-i,rAr) where r, is the row in which the i-th number appears in 
the Young tabloid. Based on Yamanouchi symbol we use the notation tn{ri,r 2 , ■ ■ ■ ,rAr_i,rAr) 
to uniquely represent a Young tabloid of partition n. 

As an example, the Young tabloid depicted in figure for partition (2,2) represents the 
equivalence class containing the four tableaux presented in figure 

Let n = [ni,n 2 ,...] and n = [n^jn^,...] be two given partitions of N (i.e. n h iV and 
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n h N) then n dominates n if for all z > 1, the sum of i greatest parts of n is greater than or 
equal to the sum of i greatest parts of n . In other words, 

i i 

n>n if and only if E Uj > E Uj for all z > 1. 

i=i i=i 

Note that in above definition of partition dominance, partitions n and n are extended by 
additional zero parts at the end as necessary. The dominance relation between two sequences 
of numbers is also known as majorization m- In terms of the Young diagrams, the number 
of squares in the first i rows of the Young diagram of partition n is greater or equal to that 
of partition n . The diagram for dominance relations between partitions of a given number is 
known as the Hasse diagram, and it is used to represent partially ordered sets. As an example 
for N = 6, partition (3,3) dominates partition (2,2,1,1) but partitions (3, 3) and (4, 1, 1) are 
incomparable, since neither dominates the other. The Hasse diagram for all possible partitions 
of Y = 6 is depicted in figure 

For n\- N , the M"" is the vector space over real numbers M whose basis consists of a set of 
tabloids of rz-shape given by 

M" = R{{tn{ri{l),r2{l),-- • YAf(l))},-'' , {'tniri{u),r2{n), ■ ■ ■ ,rAr(zz))}} 

where the set {{tn(n(l)Y 2 (l), •• • YAf(l))},-- - ,{tn{ri{i'),r 2 {n), ■ ■ ■ ,rAr(i/))}} is a complete list 
of distinct tabloids of n-shape. The symmetric group Sjy acts transitively over this set, i.e. by 
permutation the Yamanouchi symbols of a tabloid can be transformed to those of any other 
tabloid from the same set. Thus is a representation of Sjsr called the permutation module 
corresponding to n. 

For n = (ni,rz 2 ,--- ^uk) k Y, the Young subgroup of Sn corresponding to n is defined 
as Sn = 5'ni X Sn 2 X • • • X SriKi where permutes 1, 2, • • • ,ni, Sn 2 permutes rzi + l,ni + 
2, • • • , m + 712 and so on. The order of the Young subgroup of n-shape is ni!n 2 ! • • • uk- Since Sn 
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( 6 ) 

t 

(5,1) 


t 

(4,2) 



( 2 , 2 , 1 , 1 ) 


\ 

( 2 , 1 , 1 , 1 , 1 ) 

i 

( 1 , 1 , 1 , 1 , 1 , 1 ) 

Figure 4: The Hasse diagram for all possible partitions of T = 6. 

is a subgroup of the number of left or right cosets of Sn 'm Sn is N\/{nil • n 2 ! • • • nxO which 
is also number of distinct tabloids of n-shape or dim{M"‘), hence there is a bijection between 
and the {IljP}, where {Hj} is a transversal for Sn in Sn- 

2.3 Cayley Graph &; Schreier Coset Graph 

Let ^ be a group and let 5 T TL. The Cayley graph of Tl generated by S (referred to as the 
generator set S), denoted by Cay{Tl,S), is the directed graph G = (V,T) where V = TL and 
£ = {(x,xs)|x £TI,s £ 5}. If 5 = 5“^ (i.e., S is closed under inverse), then Cay(Tl,S) is an 
undirected graph. If TL acts transitively on a finite set 11, we may form a graph with vertex set 
V = n and edge set £ = {{n, ns)\i' £ s £ 5}. Similarly, if Q is a subgroup in TL, we may form 
a graph whose vertices are the right cosets of Q , denoted {TL : Q) and whose edges are of the 
form £ = {(Qh, Qhs)lQh £ (TL : Q), s £ S}. These two graphs are the same when 17 is the coset 
space (TL Q), or when Q is the stabilizer of a point of 17 and is called the Schreier coset graph 
Sch{n,S,Q). 
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3 Classical Continuous Time Consensus (CTC) 

Consider a group of N agents with an underlying connected graph Q = (V,£l). Each edge {i,j} 
indicate bidirectional communication between agent i and agent j , thus the resultant underlying 
graph Q is an undirected graph. 

Let Xi be the state of agent i. In the continuous time consensus (CTC) algorithm, each 
agent’s dynamics evolves according to the following state dynamics equation, 

N 

Xi{t) = '^Wij{xi{t) - Xj{t)), for i = l,...,N, (2) 

Defining the vector x = [xi, • • • ,X]^Y the vector of states, we can rewrite the above state 
evolution formula in compact notation as below, 

x{t) = —Lg X x(t), t G M > 0, (3) 

where Lg is the graph Laplacian matrix for a weighted graph as defined in section |2.1[ It is 
well-known that limt_>.oo —)> 1 x I'^/N where 1 is the left eigenvector of Lg corresponding 

to eigenvalue 0. Thus according to the state dynamics equation Q it can be concluded that 
limt^oo x{t) —)• X X x(0) = ^ Xj(0)l. In other words, the final equilibrium state of the 
consensus algorithm is the average of agents’ initial states if the underlying graph is connected. 
The convergence rate of the algorithm to its equilibrium state is governed by the second smallest 
eigenvalue of the graph Laplacian (A 2 ) [23]. Larger values of A 2 results in faster convergence rate. 
The CTC algorithm is also known as the continuous time Markov chain algorithm. The inverse 
of the second smallest eigenvalue of the graph Laplacian (A 2 ) is referred to as the relaxation 
time [T7] . 

For a given connected network with an underlying graph topology the Fastest Continuous 
Time Consensus (FCTC) problem can be formulated as below, 

max \2{Lg) 

W 

^ ( 4 ) 

s.t. 2_^ Wjk < D, 

{j,k}££ 

where Wjk is the weight on the edge from node j to node k and D is an upper limit on the 
total amount of weights. This optimization problem can be defined in the form of standard 
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semidefinite programming (see Appendix [B|) as below 

min — s 

W 

S.t. Lg + ( —s)/ — (1 X l'^)/N ^ 0, 

D - ^ Wjk > 0. 

{j,k}&£ 

In the formulation above 1 is the column vector of all one. 

To the best of our knowledge, analytical optimization of the CTC problem has been addressed 
only for tree topologies in m by algebraic method. In section we have provided analytical 
solution to the semidefinite programming formulation of the CTC problem for a wider range of 
topologies. 

An automorphism of the graph Q = (V,T) is a permutation cr of V such that {i,j} G T if 
and only if {a{i),a{j)} G T, the set of all such permutations, with composition as the group 
operation, is called the automorphism group of the graph and denoted by Aut{Q). For a vertex 
i G V, the set of all images cr(i), as a varies through a subgroup G C Aut{Q), is called the orbit 
of i under the action of G. The vertex set V can be written as disjoint union of distinct orbits. 
In [12] , it has been shown that the optimal weights on the edges within an orbit are equal. 

4 Continuous Time Quantum Consensus 

4.1 Lindblad Master Equation 

We consider a quantum network as a composite (or multipartite) quantum system with N 
qudits. Assuming TL as the d-dimensional Hilbert space over C, then the state space of the 
quantum network is within the Hilbert space ® The state of the quantum 

system is described by its density matrix (p). This matrix is positive Hermitian and its trace 
is one {tr{p) = 1). The network is associated with an underlying graph Q = {V,T}, where 
V = {1,...,A^} is the set of indices for the N qudits, and each element in £ is an unordered 
pair of two distinct qudits, denoted as {j, A;} G T with j, /c G V. Permutation group Sn acts in a 
natural way on V by mapping V onto itself. For each permutation tt G we associate unitary 
operator U-^ over TL®^, as below 

Cr(^l ® * * * ^ Qn') 1 ^ 77 ( 1 ) ^ * * * ^ Q'k{N)^ 
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where Qi is an operator in % for all i = 1,..., A special case of permutations is the 
swapping permutation or transposition where 7r(j) = A:, TT{k) = j and 7r(i) = i for all i G V and 
i ^ j,k We denote the swapping permutation between the qudits indices j and k by -Kj^k and the 
corresponding swapping operator by Uj^k- la Appendix [A| the swapping operator Uj^k has been 
expressed as linear combination of the Cartezian product of Gell-Mann matrices. 

Employing the quantum gossip interaction introduced in [21], the evolution of the quantum 
network can be described by the following master equation 


dp 

dt 




X p X - p 


( 5 ) 


where Wjk is the positive constant weight over the edge j, k. These weights form the distribution 
of limited amount of weight up to D, among edges of the underlying graph, i.e. 


^hk < D. ( 6 ) 

In order to have the set of transpositions corresponding to the edges of the underlying graph as 

the generator set S of the symmetric group Sn, the underlying graph should be connected. 

For the case oi H = 0, the evolution of p(t) is described in the following Lindblad master 
equation 

'^j,k(UjkX pxU^.^^- py ( 7 ) 

which is named Quantum Consensus Master Equation (QCME) by authors in |27|, and its 
resultant quantum consensus state m is defined as 

'’* = M E u.pmi (8) 

In |27] . it is shown that the QMCE reaches quantum consensus, namely limi^oo p(t) = P* 
provided that the underlying graph of the quantum network is connected. 

The aim of the analysis presented in the rest of this paper is to evaluate and optimize the 
convergence rate of the QCME to its quantum consensus state. To this aim, we expand the 
density matrix (p) as the linear combination of the generalized Gell-Mann matrices (introduced 
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in Appendix [A]) as below, 


1 


d2-l 


^ 2^ ^ y Pui,U2,---,Un ■ ® 

Ui,U 2 ,-",Un =0 


'Un ; 


(9) 


where N is the number of particles and (8> denotes the Cartesian product and A matrices are 


the generalized Gell-Mann matrices as in (74) and (76). Note that due to Hermity of density 


matrix, its coefficients of expansion numbers and because of unit trace of p 

we have po,o,...,o = 1- 

Using the decomposition of p in ([^, its permutations can be written as below 


d^-i 


Uj^k X p X ® ® 


X 


UN 


( 10 ) 


Ui,U2,---Un=0 


Note that in (10) due to permutation operators, the place of indices pj and pk in the index of 
parameter p are interchanged. Substituting the density matrix p from (© and its permutation 
(10) in Lindblad master equation ([^ and considering the independence of the matrices 
A,, 


'Ml 


'M 2 


• • A^jy we can conclude the following for Lindblad master equation ([^ , 
d 


dt 


7*M1 i--- iMjv 


{Pui,— ,Uk,—,Uj,—,U n Pui,— ,Uj,—,Uk,—,U n) 

{j,k}££ 


( 11 ) 


for all pi, p2, • • • , PAf = 0, • • • , - 1, 


with the constraint ^ on the edge weights. Following the same procedure, the tensor component 
of the quantum consensus state (j^ can be written as below 


P 


* 


1 

m 


^ir(Ml).n'(M 2 )v,ii(Miv)(l*) 

ttSSjv 


( 12 ) 


and for the connected underlying graph, the QCME reaches quantum consensus, componentwise 
as below 

}}:^Pui,U2,---,UNi^) — Pui,U2,---,UN^ 

Comparing the set of equations in (0 with those of the CTC problem in (|^ we can see that 
the Quantum Consensus Master Equation 0 is transformed into the classical CTC problem 0 
with d"^^ — 1 tensor component P/xi,---,/ijv as the agents’ states. Defining Xq as a column vector 
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of length with components the state update equation of the classical CTC can be 

written as below, 

( 13 ) 

Lq is the corresponding Laplacian matrix as below, 

Lq= ^ Wj^k{I^2N - Uj^k), ( 14 ) 

{j,k}&£ 

where Uj^k is the swapping operator given in Appendix [A| (|7^ , provided that d is replaced with 
which in turn results in Gell-Mann matrices of size d‘^ x d?. As explained in section the 
convergence rate of the obtained CTC problem is dictated by the spectral gap of its associated 
underlying graph which is the second largest eigenvalue {X2{Lq)) of its Laplacian matrix Lq. 
Thus, the corresponding Fastest Continuous Time Consensus problem can be written as the 
following optimization problem. 


max \ 2 {Lq) 

W 

s.t. Wjk < D. 

{j,k}£S 


(15) 


We refer to this problem as the Fastest Continuous Time Quantum Consensus (FCTQC) prob¬ 
lem. This is the same optimization problem as in 127! with the difference that in |27| it has been 
obtained in computational basis. 

The QCME Q reaches quantum consensus ([^, due to the fact that the generating set is 
selected in a way that the whole group of S'at can be generated, and the resultant Cayley graph 
of S]\f is connected. Even though, the quantum consensus is achieved but surprisingly, the 


equations in (11) indicate that all agents are not able to exchange information with each other. 
This is due to the fact that the underlying graph of the CTC problem obtained from 0 is not 
connected and the consensus is not reachable in the same sense as in the classical CTC problem, 
where the sufficient condition for reaching consensus is the connected underlying communication 
graph. 

From the right hand side of the equation ( [IT| ), we can see that the tensor components 
that can be transformed into each other by permuting their indices are communicating and 
exchanging information with each other. These tensor components or agent states correspond 
to Young tabloids of the same partition that in turn are equivalent to the agents in the classical 
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CTC problem. Thus, the agents belonging to the same partition form the connected components 


of the underlying graph of the classical CTC problem (11). 


As mentioned above the underlying graph is a cluster of connected components where each 
connected graph component corresponds to a given partition of N into K integers, namely 
N = ni + n 2 + • • • + uki where K < (P and nj for j = 1,... ,K is the number of indices in 
Pm,U 2 ,--;UN with equal values. For a given partition and its associated Young Tabloids, more 
than one connected component can be obtained depending on the value of the fi indices. As an 
example consider a quantum network with three qubits and the path graph as its underlying 
graph. In this network, the values that the fj, indices can take are 0,1,2 and 3. For partition 
n = (2,1) and Young Tabloids tn,(l, 1, 2), t„(l, 2,1), t„(2,1,1) and /ii = 0 and /r 2 = 1 the 
obtained underlying graph of the CTC problem is a path graph with three vertices where each 
vertex corresponds to one of the Young Tabloids mentioned above. Now for the same partition 
and Young Tabloids but different values of the ju indices (e.g. = 1 and /X 2 = 0) the obtained 

underlying graph of the CTC problem is same as that of the previous example. As a matter 
of fact for this partition, there are 12 connected components which are identical to each other. 
Each one of these connected components has an identical impact on the convergence rate of the 
QCME to its quantum consensus state. Therefore for each partition we consider only one of 
them, and we refer to this graph as the induced graph. The only exception is the case of A^ = 
where there is only one connected component corresponding to the partition that all indices are 
different from each other. These induced graphs are the same as those noted in m- 

Eor the given partition n = [ni,n 2 , • • • jUk], using the Yamanouchi symbol (introduced in 


section 2.2) a Young tabloid of partition n is uniquely represented by the notation tn{ri,r 2 , ■ ■ ■ , tn-: 


Each Young tabloid tn{ri,r 2 , • • • , r^) is equivalent to an agent in the induced graph of the CTC 
problem and its corresponding coefficient ,... ) is equivalent to the state of that agent. 

The CTC equation obtained from 0 for partition n is as below, 


(m) V iP'T'jY (m) 


(16) 




where m varies from 1 to = A^!/(ni! • 712 ! • • -nK^-) and TTj^i transposes the j-th and l-th Ya¬ 
manouchi symbols i.e. rj and ri. Note that for the agent states that their Yamanouchi symbols 


TAf)- 
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{rj,ri) are equal, the value inside the summation above is zero. 


We define the column vector as the state vector of the associated CTC problem (20) 


of agiven partition n. This vector includes the tensor components corresponding to the Young 
Tabloids of the partition n and it has iz = N\/(nil ■ n 2 C ■ ■ nx !) elements. As mentioned above the 
underlying graph of the CTC problem is a cluster of connected components, i.e. the Laplacian 
matrix Lq is a block diagonal matrix where each block corresponds to one of the connected 


components, with state vector Xn- The state update equation (13) for the state vector Xn is 
as below, 


dXr 

dt 


— Ln^ni 


(17) 


with Ln as the Laplacian matrix which is one of the blocks in Lq. 


The tensor component of the quantum consensus state (12) for partition n takes the following 
form 


N\ E 


jY! ^M7r(ri)v?M7r(r;y) 

7rG5iv 


(18) 

As explained in section [2^ Sx acts transitively over the set of Young tabloids or agents and 
consequently over the following set of agent states 


{{{pn{r i{l), r 2 {l)x ■ ■ ,r7 v(l))},-- - Apn{ri{v),r 2 {v),-■ ■ ,rw(i^))}}) 


with the Young subgroup Sn as its stabilizer subgroup. Since the group elements of the Young 
subgroup do not change the Yamanouchi symbols. Based on the one to one correspondence 
between agent states and the right or left cosets of Sn in Sx, it can be concluded that the 
connected component is the Schreier coset graph of permutation group Sx with Young subgroup 
Sn and generating set consisting of transpositions associated with edges of the underlying graph 
of the quantum network. For the case of trivial Sn (i.e. n = [1,1,...,!]) the Schreier coset 
graph is reduced to Cayley graph. 

In the following, we provide some of the typical partitions with their corresponding connected 
graph component. 

The most simple case is the one that all indices are the same, i.e. the partition is n = (ni) = 
{N). The Yamanouchi symbols for this partition are ri = r 2 = • • • = rw = 1 and the CTC 
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equation (16) is as below 


d 

71 


= 0 , 


for = 0,1 ,..., — 1. 


and therefore 


= dw,w,-,w(0), for iai = 0,l,...,d^ -1. 


The induced graph of this partition is the edgeless or the empty graph that is a graph without 
any edges. This is the Schreier coset graph Sch{SN,<S, S^). Due to lack of any information 
exchange between agents, the agent states does not change by time, and they maintain their 


initial values. Thus for the quantum consensus state (18) we have 




IV! E d7r(/ii),7r(/ii),...,7r(/ii)(0) PfJ-l (0), 


7rg5jv 


where the second equality above is based on the fact that the agent state remains 

intact under the permutation vr or any exchange of information. 

The next non-trivial partition is the case where all p indices are the same except for one of 
them, thus we have 

n = [iV- 1,1], 


i.e. the Yamanouchi symbols are as = 1 for i = {1,..., N} \ {j} and rj = 2. Thus the agent 
state can be written as where for ease of notation we denote the agent state 


by the scalar variable Xj for j = Hence the CTC equation (16) for the partition 

n = [A^ — 1,1] can be written as below. 


j^Xj= Wj^k{xk-xj), ( 19 ) 

with the constraint ([^ on the edge weights. M{j) is the set of neighbours of node j in the graph 
Q. Considering Xj as the state for node j, the equation above is same as the classical CTC 
problem over the underlying graph Q which in turn is the Schreier graph Sch{Si\f,S, Sjsf-i). For 
this particular partition, the induced graph of the partition is same as the underlying graph of 
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(a) 


(b) 


Figure 5: The induced graphs for a path graph with 3 vertices for partitions (a) n = (2,1) and 
(b) n = (1,1,1). 


the quantum network {Q). For the quantum consensus state (12) of this partition we have 


1 iV ^ N 

M2 ,M1v,Mi(Q) = Ti (Q) i 


N 


i=i 


2-th 


i=i 


Note that, in this case, the quantum consensus state is same as the final equilibrium state of 
the classical CTC problem. 

For the case that all indices are different, namely for the partition n = [1,the 
CTC problem is referred to as interchange Process [6]. This case is possible if < (P. The 
Yamanouchi symbols for this partition take different values from 1 to where no two symbols 


are equal to each other. The quantum consensus state (12) for this partition is same as (18) with 


the exception that no two p, indices have the same value. The induced graph of this partition 
is the Schreier coset graph Sch{S]\f,S,e) where e is the identity element of Sn- This Schreier 
coset graph is same as the Cayley graph (Sat, 5). 

As an example, consider the path graph with three vertices (denoted by Gps) as the un¬ 
derlying graph of the quantum network. For partition n = (2,1) over graph Qp^ the induced 
graph is as depicted in figure]^ (a) which is same as the underlying graph Gps- But for partition 
n = (1,1,1) the induced graph obtained is a cycle graph as depicted in figure]^ (b). 
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4.2 Intertwining of Induced Graph 


As stated before, the aim of this paper is to evaluate and optimize the convergence rate of the 
QCME Q to its quantum consensus state Q. In doing so in section 4.1 we have shown that the 
QCME Q can be modeled as a classical CTC problem with a disconnected underlying graph 
that is a cluster of connected components referred to as the indnced graphs. Therefore, the 
convergence rate of the QCME Q is equivalent to the convergence rate of the obtained classical 
CTC problem, i.e. the second largest eigenvalue of its Laplacian matrix (A 2 (Tq)), which in 
tnrn is determined by the second largest eigenvalues of the induced graphs. In this section, we 
address the order of the second largest eigenvalues of the induced graphs. We have shown that 
the second eigenvalues of all induced graphs are equal and based on this we will show that the 
general problem reduces to finding the second largest eigenvalne of the induced graph of the 


least dominant partition in the Hasse diagram (as explained in section 2.2). To this aim, we 


show that the eigenvalues of the induced graph corresponding to the dominant partition (higher 
level of the Hasse diagram) is included in the eigenvalnes of the less dominant partition (lower 
level of the Hasse diagram). 

For ease of notation we consider the Young Tabloid where the Yamanouchi symbols are 
sorted i.e. = 'r £)._|_2 = • • • = roi+m = i with Di = = 0- is obvious 

that the other Tabloids of this partition can be obtained from the permutation of Yamanouchi 


symbols of the above Young Tabloid. Using this notation, the CTC equations (16) for a given 
partition can be written as below 


A 

dt 


Pu 7r(l) v)/^7r(l) v!p'7r(/C) v?/^7r(iY) 




;7r{l) i-'-iMiTj- J7r(l) viA'tTj' ;ir(K)v,A‘7rj_;7r(K) 


( 20 ) 


Pu 


where vr G permutes the location of indices and transposes the location of j-th and Z-th 
indices. 

In Hasse diagram, the one level dominance i.e. the partitions that are one level apart in 
Hasse diagram can be classihed into two categories. First category is the case when one of the 
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boxes in the Young diagram is displaced from a higher row to an existing lower one in the Young 
diagram, provided that the new diagram is again a Young diagram. An example for this category 
is the partitions (4, 2) and (3, 3) in the Hasse diagram depicted in figure]^ The second category 
is the case when one of the indices is changed to a new value. In terms of Young diagram of 
the partitions, one box is moved from a higher row to a new row at the bottom of the Young 
diagram. An example for this category is the partitions (4, 2) and (4,1,1) in the Hasse diagram 
depicted in figure 


First Category 


Consider two given partitions of N namely, n and n , where partition n is one level dominant 
to partition n . If the dominance level is of first category then Partition n can be written in 
terms of partition n = [ni, n 2 ,..., n^] as following n = [ni, n 2 ,..., Um — 1, , nr-i,nr + 

1,... jUk], provided that Um > ?^m+l and < n^-i. Considering the derivation of the CTC 


equation (20) for partition n , we can define the new variable p in terms of tensor components 
p as below 


P Un(l) (1) 5”'5/^7r(m)v!/^7r(m)5*--)/^7r(r)v5M7r(r)i’”)/^7r(it)!"*5/^T 


(K) 


Pu 7r(l) • * )/^7r (l) • ■ i/^7r(r) (m) v ■ 5/^7r{m) v * (r) • • jMtt( r) ) ■ ■ • 5/^7r( JY) v * • (iY) 

S---^ S---^ V---^ V---^ 


+ Pu 


7r(1) r • • >/^7r(l) r • • ?/^7r(m) ?/^7r(r) i/^7r(m) ? • • • ’/^tt( m) v • • )f^7r{r) v ?/^7r(r) >/^7r(r) r • • ?/^7r(iY) v • if^7r{K) 


( 21 ) 


+ 


+ Pu 


7r (1) r • • >/^7r(l) r • • ?/^7r(m) v »/^7r(m) !M7r(r) 5 • • ■ ?/^7r(r) ! • • • !M7r(r) r • • ?/^7r(iY) r • • ?M7r(iY) 


Taking the derivative of p in (21) and applying the CTC equation of partition n (20) to the 


right hand side of the resultant equation, it is straightforward to show that p obeys the same 


CTC equations of partition n as in (20) 


Second Category 


If the dominance level is of second category then Partition n can be written in terms of partition 
n = [ni, 712, • • •, nK] as = [ni,..., Um — 1, Um+i, • • ■, Ufc, 1], provided that Um > ?^m+i- Similar 


to first category, considering the derivation of the CTC equation (20) for partition n , we can 
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define the new variable p in terms of tensor components p as below 


Pu 7r(l) ! • • • iM7r (1) V >M7r(m) v ?/^7r(Tn) 5 • • • 


Pu 7r(l) !■ ■ • i/^TT(1) )• *■ )M7r(m) v• 5/^7r(m) )• *• (i^) )■ ■ ■ i/^7r( ) Mtt( iil + 1) 

V > V ^ V w > 

^ v ^ v ^ V ^ v' 

nr<' 1 


Tim —1 


+ Pm. 


m) v?M7r(m) iM7r(i<'+l) v?M7r(i<') i'“iI^'k{K) /^7r(m) 


Tim-l 


+ Pm 


m) V? M7r(m) ’Mtt( i^’+l) jMtt( m) r • • r • 't^Tz^K) /^7r(m) 


( 22 ) 


+ 


+ Pm 


■7r(l) ’•••?M7r(l) viM7r(X+l) i^Tz{m) 1•••1^^z{m) ^’••^^7z(K) ^‘■■>l-^7z{K) ^7z(m) * 


Tim-1 


In the same manner as in first category, after taking the derivative of p in (22) and applying 


the CTC equation of partition n (20) to the right hand side of the resultant equation, it is 


straightforward to show that p obeys the same CTC equations of partition n as in (20) 


For both categories, we have shown that the newly defined variable p obeys the same CTC 


equations of partition n as in (20). Using this result, in the following we show that all eigenvalues 


of the partition n are amongst the eigenvalues of the partition n . 

we define the column vector Xn as the state vector for partition n with p as it compo¬ 


nents. Depending on the category p can be from the left hand side of either equations (21) or 


(22). Recalling X / as the state vector of partition n defined in section 


4.1 


in terms of tensor 


components p, we can conclude the following, based on equations (21) and (22) 


Xn = p{n n) X X 


(23) 


where P{n —)• n ) is the projection matrix with u = N]/(nil - 712 ! • • • n^O rows and u = N\/(ni\ ■ 


Un! 


I.... 


2 ; • • ■Ti^/!) columns. Matrix P(n —)• n ) is the projection matrix for the surjective projection 


that maps the states of partition n onto sates of partition re. Taking the derivative of the 


equation (23) and substituting the derivative of the state vectors with the products of state 


vectors and their associated Laplacian matrices according to state update equation ( [l7| ) we 
can conclude that the following relation holds between the Laplacian matrices associated with 
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partitions n and n 


Ln X P{n —)> n) = P{n —)• n) x L^/, (24) 

This is known as the intertwining relation. By taking the transpose of both sides of this equation 
we obtain 


P^{n —)• n) X Ln = L^i x P^(n —)■ n), 

Then for each eigenvalue of Ln denoted by 7 and its associated eigenvector T we have 


(25) 


X r = 7 • r. 

Multiplying both sides by P^(n —)• n ) from left 


(n —)> re ) X X r = 7 • I P^ (re — >■ re ) x T I , 


Using (25) we obtain the following 


X (^P'^in ^ re') X r) = 7 • (^P^(re ^ re') x t) , (26) 

Since P(re —)• re') is the projection matrix for a surjective projection then P^(re —)• re') is the 
projection matrix for an injective projection and it does not have a null space i.e. it doesn’t 
have zero eigenvalue m Therefore,s it can be concluded that any eigenvalue of Ln is also an 
eigenvalue of L^/. 

Based on this conclusion and the fact that the first eigenvalue (Ai) of both Laplacian matrices 
Ln and L^/ are zero then we can conclude the following relation for the second eigenvalue of 
the Laplacian matrices, 

\2{L^>) < X2{Ln) (27) 

In simple words, we have shown that the second eigenvalue of a partition is less than or equal to 
that of its one level dominant partition. Applying this conclusion to the partitions in the Basse 
diagram of a given N, we can conclude that the second largest eigenvalue of the partition at 
bottom (top) of the Basse diagram is the smallest (greatest) and the second largest eigenvalue of 
all other partition are in between with the order same as the dominance order in Basse diagram. 
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In other words, 


[i,^) < A2([2, 1^) < • • • < A2([iV - 1,1]) 

N N-2 


(28) 


Note that the partition n = [A^] has only on eigenvalne that is zero. 

In the prominent work [6] authors have proved that the second eigenvalues (i.e. the spectral 

gap) of the partitions and [N — 1,1] (known as the interchange and the random walk 

N 

processes, respectively) are equal. This is known as the Aldous’ conjecture [2]. Considering this 


result and the relation (28) it can be concluded that the second eigenvalues of all partitions 


(except [N]) in the Hasse diagram are equal to each other. This is the generalization of the 
Aldous’ conjecture to all partitions (except [N]) in the Hasse diagram of N. 

An interesting point is a case where one of the induced graphs namely G' with N' vertices 
serves as the underlying graph for a quantum network. Since the induced graph corresponding 
to partition (N' — 1,1) is same as the underlying graph G' , therefore, the spectral gap and the 
convergence rate of all induced graphs corresponding to partitions of N' are same as those of 
all induced graphs corresponding to partitions of N. 

In the following, we provide examples of the projection matrices P{n —>■ n) for both cate¬ 
gories of one level dominance as mentioned above. For the first category, we consider the path 
graph with four vertices as the underlying graph of the quantum network. The projection matrix 
from partition n = [3,1] to partition n = [2, 2], and their Laplacian matrices are as below. 


Pin —>■ n ) = 


110 10 0 
10 10 10 
0 110 0 1 
0 0 0 1 1 1 


(29) 


Ln — 


Wl,2 -Wl,2 

-Wl,2 Wl,2 + 

0 —rC2,3 

0 0 


0 0 

-^ 2,3 0 

W2,^ + rC3,4 -t(^3,4 

-W'3,4 ^i'3,4 


(30) 
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P, 




p^ 




P. 




W 


1,2 


w. 


2,3 

(a) 


W- 


3,4 


Pnujui,P2,jui 



(b) 


Figure 6 : The induced graphs for a path graph with 4 vertices for partitions (a) n = (3,1) and 
(b) n = ( 2 , 2 ). 


W2,3 -W2,3 0 

-102,3 Wl,2 + W2,3 + W3,4 -Wi,2 

0 -Wl,2 Wi,2+W3,4, 

0 —W3,4 0 

0 0 —W3,i 

0 0 0 


0 0 0 

— W3,4 0 0 

0 —W3,4 0 

Wl,2+li'3.4 -li'1,2 0 

-Wi,2 Wi,2+W2,3+W3,4 -^ 2,3 

0 -W2,3 W2,3 


(31) 


The induced graphs of these partitions are depicted in figure 

As an example for the second category, consider the path graph with three vertices as the 
underlying graph of the quantum network. The projection matrix from partition n = [2,1] to 
partition n = [1,1,1], and their Laplacian matrices are as below, 


Pin —>■ n ) 


1 1 0 0 0 0 

0 0 1 1 0 0 , 

0 0 0 0 1 1 


(32) 
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Wl,2 


0 





L-n 

= 

-Wl,2 

^^^l,2 + W2,3 

-W2,3 

1 


(33) 




0 

-W2,3 

W2,3 




Wl,2 + W2,3 

— W2,3 


-Wl,2 

0 

0 

0 


-W2,3 

+ W2,3 


0 

-Wl,2 

0 

0 


-Wi^2 

0 

^1^2 + W2,3 

0 

-W2,3 

0 

(34) 

0 

-Wi,2 


0 

^^^l,2 + 'W^2,3 

0 

-W2,3 


0 

0 


—W2,3 

0 

Wl^2 + ^^^2,3 

-Wl^2 


0 

0 


0 

-W2,3 

-■^1,2 

^,2 + ^i'2,3 



The induced graphs of these partitions are depicted in figure 


5 Optimization of the Convergence Rate 


In the previous section we have modeled the Continuous Time Quantum Consensus problem as 
the classical CTC problem where its underlying graph is a cluster of induced graphs. Further¬ 
more, we have shown that the second eigenvalue of the Laplacian matrices of all the induced 
graphs are the same. This eigenvalue effects the convergence rate of the CTC problem and 
therefore that of the Continuous Time Quantum Consensus problem. In this section we address 
the optimization of the second eigenvalue of the Laplacian matrices of the induced graphs in the 
modeled CTC problem. In previous section, this problem is introduced as the Fastest Continu¬ 


ous Time Quantum Consensus (FCTQC) problem (15). Considering the fact that the value of 
this eigenvalue is the same for all induced graphs then this optimization problem is reduced to 
optimizing the second eigenvalue the Laplacian matrix of the induced graph corresponding to 
partition n = [iV — 1,1]. 

In the following, first we provide the optimal results for all possible connected topologies 
with N = 2,3 and 4 vertices serving as the underlying graph of the FCTQC problem. Then 
we study the FCTQC problem over quantum networks with different topologies in their general 
form. We categorize these topologies into two groups. First group, is the group of topologies 
that the FCTQC problem have been solved using linear programming [3j and the second are 
those topologies that the FCTQC problem have been solved using semidefinite programming 
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Figure 7: All possible connected underlying topologies with Ai = 4 vertices which are non¬ 
isomorphic. 

[3]. For the complete-cored star topology we have provided the detailed solution while for the 
rest of the topologies we only report the final optimal weights. 

5.1 The Optimal Results for all topologies with N = 2,3 and 4 vertices 

Here we provide the optimal weights and the second smallest eigenvalue of the Laplacian matrix 
for all possible topologies with N = 2,3 and 4 vertices which are connected and non-isomorphic. 
The order of presentation of graphs here is in terms of increasing value of A 2 for a given value 
of 17. In other words, for a given D path graph has the least convergence rate while complete 
graph has the most. 

For a network with N = 2 vertices, the only connected topology is the path graph with 
2 vertices. The optimal value of the second smallest eigenvalue (A 2 ) for path topology with 2 
vertices is 2D and the optimal weight is D. In case of a network with N = 3 vertices, there 
are two connected topologies, namely, path topology with 3 vertices and the triangular topology 
which is a complete graph. For the path topology with 3 vertices, the optimal value of the second 
smallest eigenvalue (A 2 ) is D/2 and the optimal weight is D/2. For the triangular topology, the 
optimal A 2 and weight are D and D/3, respectively. 

Path graph with N = 4 depicted in figure (a): 


Wo = 2D/3, wi = 377/10, 


(35a) 


A 2 — 77/5, 


(35b) 
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Star graph with N = 4 depicted in figure (b): 


w = D/?>, 

(36a) 

A 2 = D/2), 

(36b) 


Lollipop graph with N = A depicted in figure (c): 

w-i - D 

Cycle graph with N = 4 depicted in figure]^ (d): 

w = D/4, 

A2 = d/2, 

Paw graph with N = 4 depicted in figure (e): 

wq = 0, wi = D/4, 
A2 = D/2, 


Complete graph with N = 4 depicted in figure]^ (f): 


w = D/6, 

(40a) 

X 2 = 2D/3, 

(40b) 


Note that for Paw graph the weight on the diameter {wq) is zero and the optimal weights and 
A 2 are same as those of Cycle graph. 

5.2 Linear Programming 

In this section, we provide the optimal results for a number of topologies where the FCTQC 
problem can be solved using Linear Programming [3]. 


(39a) 

(39b) 


(38a) 

(38b) 


(37a) 

(37b) 
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5.2.1 Cartesian Product of Edge Transitive Graphs 

This topology is obtained from Cartesian product of m edge-transitive weighted graphs. The 
weighted Laplacian matrix for the whole graph can be written as below, 

pW 0 0 ... 0 ® ® ® ■■■ ® iNm 

where Lf is the weighted Laplacian matrix of i-th graph. Ijsi. is the identity matrix with size 
Nj and Nj is the number of vertices in the j-th graph. Due to edge-transitivity, all edges of 
each edge-transitive graph have the same optimal weight. Thus the weighted Laplacian matrix 
{Lf) for each one of the graphs can be written as {Lf) = Wi ■ Li in terms of its unweighted 
Laplacian matrix {Li). Using this relation the weighted Laplacian matrix for the whole graph 
can be derived as below, 

L"' = Wi ■ Ini ® In2 (8) • • • (8) lNi_i 'Si Li® iNi+i ® ■ ■ ■ ® In^ (^1) 

We denote the eigenvalues of the Uth unweighted Laplacian matrix Lj in their sorted order by 
where ai varies from 1 to W- Using this notation the eigenvalues of the weighted Laplacian 
of the whole graph can be written as below, 




— Wl • Al^cfj -|- W2 ■ ^2,a2 T ■ ■ ■ Wm ' ^m,c 


(42) 


where ai for i = 1,... ,m varies from 1 to W- Based on the derivation in (42) and considering 


the fact that hrst eigenvalue of each unweighted Laplacian matrix Lj is zero (i.e. = 0 for 

i = l,...,m), the second smallest eigenvalue of the weighted Laplacian of the whole graph 
can be written as A^ = minrci • Ai^ 2 , ^^2 • ^ 2 , 2 , • • • , Wm ■ Am, 2 - Using this result the optimization 
problem for the FCTQC problem can be written as below. 


max s = mmwi ■ Aj 2 , 

Wl,W2,--- ,Wm i ’ 


s.t. 


Ej ■ Wj = D. 
i=i 


(43) 
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where Ej = Ej • and Ej is the number of edges in the j-th edge-transitive graph. For 

ki^j 

the optimal answer we have 


S = Wl- Ai,2 = W2 ■ A2,2 = • • • = Wm ' Xm,2- 


(44) 


From this relation we can conclude the following for the optimal value of A 2 and the weights 


A 2 = s = 


D 


^ LS 


Wj = -— for j = 


(45a) 


(45b) 


where N = 

As an example consider the Cartesian product of two complete graphs each with Ni and N 2 
vertices. Considering wi and W 2 as the weights on the edges of each one of the complete graphs, 
then for the optimal results we have 


2D 

\2 = s = -, (46a) 

2N1N2 -N1-N2 

Wl = s/Ni, W2 = S/N2, (46b) 

An obvious example for the Cartesian product of two complete graphs is the Cartesian product 
of K 2 and K 3 as depicted in figure This graph is also known as the Prism graph. The optimal 
value of A 2 is 2D jl and for the optimal weights we have wi = D/7 and W 2 = 2D 121. 

In the following we provide the optimal answer to the FCTQC problem for two well-known 
edge-transitive graphs, namely the complete graph and the cycle graph. 

5.2.2 Complete Graph 

A complete graph with N vertices is a graph where each node is connected to every other node in 
the graph. Due to the symmetry of the graph all edges have the same weight, and its Laplacian 
matrix can be written as below, 


L\y = {N — l)-w-I — w-{J — I) = N- w- I — w-J. 
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Figure 8 : The Cartesian product of two complete graphs, namely, K 2 and K^. 
The eigenvalues of the Laplacian matrix for a complete graph are as following 


Ai = 0, A 2 = A 3 = • • • = = N ■ w. 


Considering the constraint on the summation of the weights Q in the FCTC problem, we can 
conclude the following for the optimal weight on the edges of the complete graph. 


N{N -1) 


w = D ^ X 2 


2D 

N -1' 


5.2.3 Cycle Topology 

In this topology N vertices are connected in form of a cycle. A cycle graph with four vertices in 
depicted in figure (d). Cycle graph is edge transitive and thus the optimal weight on all edges 
is the same. The optimal value of A 2 and the weight are as below 

w = D/N, 

^ 2(1 — cos (27r/A'))T) 
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Figure 9: Weighted Complete-Cored Symmetric star topology with p = 5 branches of length 
g = 3. 

5.3 Semidefinite Programming 

In this section we provide the optimal results for a number of topologies where the FCTQC 
problem can be solved using semidehnite Programming [3]. For complete-cored symmetric star 
topology we provide the detailed solution while for the rest of the topologies discussed in this 
section we only provide the optimal results. 

5.3.1 Complete-Cored Symmetric Star Topology 

Complete-Cored Symmetric (CCS) star topology with parameters {p, q) consists of p path 
branches of length q, referred to as tails. Each one of the path branches contains q edges. 
Tails are connected to each other at one end to form a complete graph in the core. A CCS star 
graph with parameters p = 5, g = 3 is depicted in figure 

Automorphism of the CCS star graph is Sp permutation of tails. Using the stratification 
method it has q+1 edge orbits, namely edges connecting, vertices in the complete part to each 
other and q class of edges in tails. Thus it suffices to consider just q+1 weights wq, wi, - ■ ■ ,Wq {as 
labeled in figure . Defining the weight matrix accordingly and using the proper orthonormal 
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basis (as introduced in m) the weight matrix transforms into a block diagonal matrix where 
the diagonal blocks in the new basis are either one of the following matrices, 


Ln — 


wi —wi 0 

—rCl Wl + W2 —W2 

0 -W2 W2 + W3 


0 

0 


0 


: • . -Wq 

0 -Wq Wq 


Wl + p ■ Wo —Wl 


0 


Li = 


—Wl Wl + W2 —W2 

0 —W2 W2 + W3 


0 

0 


(48) 


: : : •. -Wq 

0 • • • 0 —Wq Wq 

Considering the relation LI = Lq + nwoeo x Bq between matrices Lq and Li and using the 
Courant-Weyl inequalities theorem i, m the following corollary for the eigenvalues of Lq and 
Li can be concluded, 


Ai(Lo) < Ai(Li) < A 2 (Lo) < • • • < Xq-i{Lo) < Ag_i(Li) < Aq(Lo) < Aq(Li). (49) 


It is obvious from above relations that the second eigenvalue of the original Laplacian matrix 
A 2 (L) is the smallest eigenvalue of Li (i.e. Ai). 

We define the following column vectors (each with g + 1 elements) as basis for matrix Li 


(48) 


eo(j) 


1 for j = 1 

< 

0 Otherwise 


— 1 for j = i 


ei{j) 


1 for j = i + 1 


for i = 1,2 ,... ,q. 


0 Otherwise 
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Thus Li can be written as below, 


g 

Li = p ■ Wo ■ Bq X Bq + 2'^Wi ■ Bj X Bj 

i=i 

Based on the results above we can express the Fastest Continuous Time Consensus (FCTC) 
problem in the form of the semidefinite programming ( as described in Appendix [B|) as below, 


min — s 

Wi 


S.t. Li — S ■ /q+l > 0, 

i=l 


(50) 


In order to formulate problem (50) in the form of standard semidefinite programming (78), we 
define Fi,B and x as below. 


where 


Fn = 


0(g+l) 0 

0 D 


Fs = 


-I{q+1) 0 

0 0 


F — 

wo — 


p-BoX bI 0 


0 


F — 

■L 7/1 ■ - 


2 • e,- X e J 


-p 


p(p-i) 


for j = l,...,g. 


X = 


Wo rci W2 . ■ • Wq s 


c(i) = < 


0 for j = 

— 1 for j = q + 2. 


F{x) 


g 

= Fo + s ■ Fg + Wo ■ F^q + Wj ■ Fyp. > 0 

i=i 
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The dual problem can be written as below, 


max —Tr[FQxZ], 

z 

s.t. Z ^0, 

Tr[Fs X Z]= Cq+2 = -1 

rr[F^„. X Z]= Cj+i = 0 for i = 0,1,..., g. 


(51) 


Since the set of basis {e^\pL = 0,1,..., g} are linearly independent one can introduce their dual 
as {e^l/i = 0,1,..., g} such that x e^ = The dual variable Z can be written in 

terms of the dual basis as below, 


Z = 


Eu=o X ep 


ZD 


EUo X ej 


ZD 


(52) 


Using the expansion (52), the dual constraints reduce to the following 


rr[F^„xZ]=0 ^ p\zo\^ = 0, 

Tr[F^^xZ]=0 2\zjf - p\zDf = 0, for j = l,2,---,q, 


(53a) 

(53b) 


where we can conclude that. 


P ■ Izol"^ = (p -'i-) ■ \zj\‘^ , for j = l,2. 


(54) 


Using dual relations (54) the complementary slackness condition (80) reduces to 


ZD \D- 


1=1 

p{p - 1 ) 


II 

M 

(55a) 

/^=0 


Q \ 


M 

II 

o 

(55b) 


1=1 


Setting = 0 is not acceptable since it will lead to the trivial case of zq = zi = ■ ■ ■ = Zq = 0, 


thus equation (55b) reduces to 


p{p — 1 ) ^ 

D - - - Wo -p2^Wj = 0 

1=1 


(56) 
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By multiplying both sides of (55a) by vectors for z/ = 0,1, • • • , g we have 


pGy^QWQZQ + 2 GyjWjZj = szu, for = 0,1,2,q, (57) 

i=i 

where G is the gram matrix and it is defined as x x = Gn^^. The gram 

matrix can be written as below, 


G = 


1 


\/( 2 ) 


\/( 2 ) 

0 


0 -k 1 


Based on (54), we can assume that zq = ■ zd, zj = zd for j = 1,2, - ■ ■ q. Substituting 


these values in (53) we obtain the following for the optimal weights, 


Wo = s- 


(p-l)(G ^)oo + Vp(p-1 )ELi [G 1) 


0,i 


Wi = S 


P{P - 1 ) 

- 1) + pELi 

2p 


for j = l,2. 


(58a) 

(58b) 


Substituting the obtained optimal weights (58) in equation (56) we obtain the following equation 


for the optimal value of the second smallest eigenvalues A 2 , 


X 2 = s = 


2D 


(59) 


(? - 1) + 2^^(?VT)I;Ll (G-\,+pj:h., (G-'),, 


where G~^ is inverse of Gram matrix. Note that matrices G and G~^ are square matrices of 
dimension q+1 but we have started the index of elements from 0 (e.g. q). This notation 

has been used due to the fact that the element in the first row and hrst column of the Gram 
matrix G refers to the inner product of the vector and cq. Thus the notation (G”^)^ ^ refers 
to the elements on first row and first column of the inverse Gram matrix (G~^). The inverse 
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Gram matrix (G is as below, 


(G 






= 


q + 1 for j = z = 0 
\/2 (g-j + l) for z = 0 , j = 1,... 
y/ 2 {q-i + l) for j = 0 , z = 1 ,..., g 
2 max(q — z + 1, g — j + 1) for z,j = l,...,g 


(60) 


Substituting the inverse Gram matrix in (59) we obtain the following results for the optimal 


value of the second smallest eigenvalue A2 and the optimal weights, 


2.0 (2p-2 + q^2p{p - 1 )) 

Wo = -7- , -V 

P{P - 1) [^P - 3 + 2q-s/2p{p - 1) + 2 pg 2 _g pgj 

2,D (^^J2p{p - l){q - j + 1) +p{q-j + l){q + j)) 

jjj ■ ^ ^_ 1 _ 

3p{q + l) (p - 1 + qy'2p{p - 1)^ + p'^q{q + l){2q + 1) 
for j = l,...,q, 


3{p - l)(g + 1) + 3^2p{p - l)q{q + 1) + pq{q + l){2q + 1) 


(61a) 

(61b) 


(61c) 


A special case of the GCS star topology is the path topology with even number of vertices 
which is obtained for p = 2. The optimal weights and the optimal value of the second smallest 
eigenvalue A2 for the path topology with 2 {q + 1) vertices are as below. 


3D(^q + 1 ) 

i2q + 3){2q + l) 

3D {(q+l)^ -f) 

; • — _ ^ ' 

' {q+l){2q + l){2q + 3) 


{q + l){2q + l){2q + 3) 


for j = l,...,q, 


(62a) 

(62b) 

(62c) 


This result is in agreement with that of Fiedler in m for a path with even number of vertices, 
(i.e by substituting p = 2 and D = 2q + l m. im), which results in 


Wo 


3(g + l) 

2q -\- 3 


3{{q + lf -f) 

j 1 

^ {q + l){2q + 3) 

6 

{q + l){2q + 3y 


for j = l,...,g, 


(63a) 

(63b) 

(63c) 
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Figure 10: The weighted graph of CCS star graph with two types of branches with p = 5 
branches of length qi =2 and q 2 = 3. 

5.3.2 CCS Star with two types of branches 

The Complete-Cored Symmetric (CCS) star with two types of branches is identified with pa¬ 
rameters {p,qi,q 2 )- This topology is a CCS star topology where two types of tails (each with 
qi and q 2 edges) are connected to each node in the complete core. A CCS star graph with two 
types of branches with parameters p = 5, qi = 2 and g'2 = 3 is depicted in figure [T0| The optimal 
weights for this topology are as below, 


,_ 92 

(p-1) {G-^)^^^ + y/p{p-l) 


i=-qi 


Wq = S X 



(64a) 


p{p - l) 



(64b) 


Wj = S X 


2p 


for j =-qi,-■ ■ ,q 2 , j^O, 
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and for the optimal value of the second smallest eigenvalue (A 2 ) we have 


X2 = s = 


2D 


92 


(p-l){G-')„„ + 2,^(?VT) £ (G-i)„.+f, £ (G-‘) 


i=-qi 

i^O 


^.*=-91 

j,i¥=0 




(65) 


Both the Gram matrix (G) and its inverse are square matrices of dimension qi + q2 + 1- 

As shown in figure [To| the index of the weights on first type of branches (with qi edges) starts 
from rc-gi. Therefore in our notation here we have used negative indexes to refer to the weights 
on these edges. 2In case of the Gram matrix and its reverse the index {—qi, —qi) refers to the 
element on first row and first column while the index ( 0 , 0 ) refers to the element on the (gi + l)-th 
row and {qi + l)-th column of the matrix. The Gram matrix (G) for the CCS Star topology 
with two types of branches is as below 


G 




-1/2 for j = i + l,z =-gi,... ,92 - 1, 

-1/2 for j = i - l,i = -91 + 1,... , 92 , * 7 ^ 0 ,1 
< -l/v/2 for {i,j} = {-1,0},{0,1},{0,-1},{1,0} 


1 for i = j,i =- 91 ,... , 92 , 
0 Otherwise 


Substituting the Gram matrix and its inverse in (65) we obtain the following formula for the 


second smallest eigenvalue (A 2 ) of the Laplacian matrix. 


A2 = s = 


6 D 

3{p - l)( 9 i + 92 + 1) + 3yj2p{p - l)Di +pD2 


( 66 ) 


with Di and D2 as below, 


= qi{qi + 1 ) + q2{q2 + 1 ), 

D2 = qi{qi + 1 )( 29 i + 1) + 92(^2 + 1)(292 + !)• 

5.3.3 Symmetric Star Topology 

A symmetric star graph with parameters p and 9 consists of p path graphs (each with 9 vertices) 
connected to one central vertex. This graph has 1+pq vertices and pq edges. A symmetric star 
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Figure 11: A symmetric star graph with p = 5 branches of length q = 3. 


graph with parameters p = 5 and q = 3 is depicted in figure 11 Due to the symmetry of the 
graph all edges with the same distance from the central vertex have the same weight, denoted 
by Wj for j = 1,..., q. The optimal value of the second smallest eigenvalue (A 2 ) of Laplacian 
matrix is 

A2 = 


6 D 


pq{q + l){2q + 1)' 
and the optimal value of the weights are as below, 


W-j = 


3D{q + j){q-j + l) 
pq{q + l){2q + 1) 


for j = l,...,q. 


5.3.4 Palm Topology 

A palm graph with parameters {p, q) consists of a path graph with q vertices connected to the 
central vertex of a star graph with p branches as shown in figure for parameters p = 5 and 
q = 4. This graph has p + q + 1 vertices and p + q edges. Due to the symmetry of the graph all 
edges connected to the central vertex of the star graph have the same weight (denoted by wq), 
except the edge connecting the path graph to the central vertex. The optimal answer varies, 
depending on the values of the parameters p and q. For 2p > q{q + 1) the optimal value of the 
second smallest eigenvalue (A 2 ) is as below. 


X 2 = s = 


6 D 

6 p + q{q + l){2q + 1) ’ 
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Figure 12: A palm graph with parameters p = 5 and g = 4. 
and the optimal weights are as following, 


Wq = s, 

(g - j + 1) ((g+ l)(2g + 1) +p{q + j)) 
2 {p + q+ 1) 


for j = I,... ,g. 


For 2p < g(g + 1) the optimal value of the second smallest eigenvalue (A 2 ) is as below, 


l2D{p + q+l) 

(g + i) (g + 2) (6 + g(g + Ap + T)) ’ 


and the optimal weights are as following. 


Wj = s- 


{q + !)(<? + 2 ) 

■^0 = ^-TAV’ 

2 (p + g + 1) 

,iQ - j + 1) jpjq + j + 2) + {q+ l)j) 
2 {p + q+l) 


for j = !,■■■ ,q- 


5.3.5 Lollipop Topology 

This topology is obtained by connecting a path graph (with g vertices) to one of the vertices in a 
complete graph with p + 1 vertices. By bridging vertex, we refer to the vertex in complete graph 
that is connected to the path graph. Considering the symmetry of the complete graph the edges 
in the complete graph can be categorized into two groups. The first group is those connecting 
the vertices in the complete graph other than the bridging vertex. We denote the weight on the 
edges of the first group by W-i. The second group is the edges connecting the bridging vertex to 
other vertices in the complete graph. We denote the weight on the edges of the second group by 
wq. The weights on the edges of the tail are denoted hy wi,W 2 , ■ ■ ■ ,Wq. The Lollipop topology 
is depicted in figure along with the weights assigned to the edges. The optimal answer varies, 
depending on the values of the parameters p and g. For g(g +1) < ■\/2p{p + 1) the optimal value 
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Figure 13: The weighted Lollipop graph with parameters p = 4 and q = 3. 
of the second smallest eigenvalue (A 2 ) is as below, 


A2 = s = 


12D{p + q+l) 
A 


where A is 


A = 6 {p- l){p + g + 1 ) 

+ {q+l) ( 6 qy' 2 p{p+l) + 6 {p + 1 ) + pq{ 2 q + 1) + qiq"^ - 1)) , 

and the optimal weights are as following, 

= ‘2 {p + g+l]{p + l) + iVMp + D) . 

s-Wo 

W -1 = -, 

p 

^ ^ 2 {p + q\-^l) + ^) + + j) + g + 1 ) > for j = l,---,q. 

For the case that q{q + 1) > ■\/2p{p + 1) the optimal value of w-i is zero and the Lollipop 
topology reduces to Palm topology where the optimal answer for this topology is provided in 
section 15.3.41 

5.3.6 Two Coupled Complete Graphs 

In this topology, two complete graphs each with A^i + N 2 and N 2 + vertices respectively, 
share N 2 vertices. In figure two coupled complete graphs with parameters Ni = 3, N 2 = 2 
and = 4 is depicted. Due to the symmetry of the complete graphs weights can be divided 
into five groups. W-2 is the weight on edges connecting the vertices on the left complete 


40 
















S. Jafarizadeh Optimizing the Continuous Time Quantum Consensus 



Figure 14; The weighted Lollipop graph with parameters p = 4 and q = 3. 

graph to each other and rc_i is the weight on the edges connecting the Ni vertices on the left 
complete graph to the N 2 vertices in the middle, wq is the weight on edges connecting the N 2 
vertices in the middle to each other. Similarly the weights wi and W2 are defined for the weights 
on the edges of the complete graph on the right-hand side of the topology. For the symmetric 
case where Ni = the optimal weights and X 2 are obtained for two following categories. If 
A'*! < N 2/2 the optimal weights and A 2 are as below 

^ ^ “ 4 N 1 N 2 + (iV2 - l)(iV2 - 2iVi) ’ 

W 2 = W -2 = 0 Wi = W-i = S/N2, Wo = {N2 - 2iVi)/iV|, 

If A^i > N 2 /2 the optimal weights and A 2 are as below 

A2 = D/{ 2 Ni), 

W2 = W-2 = 0 rci = rc_i = D/{2NiN2), wq = 0, 

From the optimal weights, it is apparent that for the last symmetric case where Ni > N 2/2 
the whole topology reduces to a 3-partite graph. For the nonsymmetric case where Ni 7 ^ N^, 
the optimal results are too long to report here. But interestingly in the nonsymmetric case if 
A^i > then the optimal value of the weight W 2 is zero. 
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6 Conclusions 

We have optimized the continuous time quantum consensus algorithm in terms of its convergence 
rate over a quantum network with N qudits. It is shown that the optimal convergence rate is 
independent of the value of d in qudits. By establishing the intertwining relation between one 
level dominant partitions in the Basse Diagram of integer N, we have shown that the spectrum 
of the induced graph corresponding to the dominant partition is included in that of the less 
dominant partition. Based on this result, the proof of the Aldous’ conjecture is extended to 
all possible induced graphs and it is shown that the problem of optimizing the convergence 
rate of the quantum consensus reduces to optimizing the second smallest eigenvalue of the 
Laplacian of the induced graph corresponding to partition {N — 1,1). By analytically addressing 
the semidefinite programming formulation of the reduced optimization problem, closed-form 
expressions for the optimal convergence rate and the optimal weights are provided. Interestingly, 
the optimal weight over some edges in certain topologies, namely Lollipop, Paw and two coupled 
complete Graphs are zero. This suggests that adding new edges does not necessarily improve 
the convergence rate. 

By extending the Aldous’ conjecture to all possible induced graphs, we have shown that 
the values of their second smallest eigenvalues are all equal. This shows that despite being 
disconnected and having different topologies, the induced graphs have the same spectral gap 
and therefore, the same asymptotic convergence rate. 


Appendix A Generalized Gell-Mann Matrices 

First we introduce the d-dimensional elementary matrices {e^\k,j = These matrices 

are d x d square matrices where the element in k-th. row and j-th. column is one, and other 
elements are zero. These matrices satisfy the commutation relation, i.e. 



- Si,iel 
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Based on the elementary matrices (e^) we define the following matrices 


e'] = e] + el 




(72a) 

(72b) 


1 <k < j < d 


In addition to off-diagonal generators, we dehne d — 1 diagonal generator matrices as below 


V 

Vr = 


r{r + 1) 


^e]-re;l\ 

i=i 


(73) 


The range of r varies from 1 to d — 1. In total d^ — 1 generators are defined in (72) and (73). 


Now we define the generalized Gell-Mann matrices (also known as A matrices) as below, 


\j-lT+2k-l = f^j-, 

\ 

Xj2_i = 


(74a) 

(74b) 

(74c) 


Note that the following relation exist between the traces of the lambda matrices defined above. 


tr{A^ X A„} = 2(5^,„ for /i, u G {0,1,..., - 1}. 


(75) 


In order to form a complete orthogonal Hermitian operator basis, we add Aq defined as below 


to the set of A matrices defined in (74), 


Ao = t/ 


(76) 


The swapping operator {Uj^k) in terms of generalized Gell-mann matrices can be written as 
below, 

Uj^k = 


d^-i 


lE 




k ,=0 


tN 


(77) 


® /d ® ® (8)Id (g) • • • (g) /d (g) (g)/d (g) • • • (g) Id + 


V 


d 


j-th 


k-th 
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As an example SU{2) generators (known as Pauli matrices) are given as below, 


Aq — h 


1 0 
0 1 


Al — ©2 — 62 + 6^ 


0 1 
1 0 


A 2 = ^2 = -*(4 - e?) 


0 -i 
i 0 


h = Vi =e\-el = 


1 0 
0 -1 


Appendix B Semidefinite Programming (SDP) 


Semidefinite Programming is a convex optimization problem that aims to minimize a linear 
function subject to a linear matrix inequality constraint [3]. It can be formulated as below, 


T 

mm c • X 

X 

\x\ 

s.t. F{x) = XjFj + Fq ^ 0, 
i=l 


(78) 


The formulation above is referred to as the primal problem. The minimization variable is the 
vector X. Vector c and matrices Fi are the problem parameters. The inequality F{x) ^ 0 
means that F{x) is a positive semi-definite matrix. 

Any vector x satisfying the constraint F{x) ^ 0 is called primal feasible point, and if x 
satisfies F{x) >- 0, it is called strictly feasible point. By convention, the minimal objective value 
■ X is called primal optimal value. Every primal problem has an associated dual problem. 
Dual problem is a maximization problem and it is formulated as below 


max — Tr[i^o x 

z 

s.t. Z ^ 0,Tr[Fi X Z] = Ci, 


(79) 


Here the optimization variable is the real symmetric (or Hermitian) positive matrix Z and the 
problem parameters are Z and Fi which are the same as in primal problem. Tr[A] means trace 
of matrix A. Any matrix Z satisfying the constraints in dual problem is called dual feasible (or 
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strictly dual feasible if ^ 0). 

The objective value of the primal feasible point is an upper bound on the minimal objective 
value cF ■ X. Similarly the objective value of the dual feasible point is an lower bound on the 
dual optimal value —Tr\FQ x Z] For a primal feasible point x and a dual feasible point Z we 
have, 

\x\ 

c?' ■ X + Tr[Fo X Z] = ^ Tr[Fi x Z]xi + Tr[Fo x Z] = Tr[F{x) x Z] > 0. 

i=l 

This proves that —Tr[FQ x Z] < cF ■ x and under relatively mild assumptions, we can have 
—Tr[F'o X Z] = • ;e. If the equality holds, one can prove the following optimality condition 

on X. 

A primal feasible x and a dual feasible Z are optimal, if and only if 

F{x) X Z = Z X F{x) = 0, (80) 

where x and Z are optimal ones. This condition is called complementary slackness condition. 

In one way or another, numerical methods for solving SDP problems always exploit the 
inequality —Tr[FQ x Z] < —Tr[FQ x Z] < ■ x < ■ x, where —Tr[FQ x Z] and ■ x 

are the objective values for any dual feasible point and primal feasible point, respectively. The 
difference ■ x + Tr[Fo x Z] = Tr[F{x) x Z] > 0 is called the duality gap. If the equality 
cF" ■ X = —Tr[FQ X Z] holds, i.e. the optimal duality gap is zero, then it is said that strong 
duality holds. 
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